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Porcine epidemic diarrhea virus (PEDV) has been detected 
sporadically in Italy since the 1990s. We report the phyloge¬ 
netic relationship of swine enteric coronaviruses collected in 
Italy during 2007-2014 and identify a drastic shift in PEDV 
strain variability and a new swine enteric coronavirus gener¬ 
ated by recombination of transmissible gastroenteritis virus 
and PEDV. 


orcine epidemic diarrhea virus (PEDV) and Trans¬ 
missible gastroenteritis virus (TGEV) (family Coro- 
naviridae, genus Alphacoronavirus ) are enveloped viruses 
that contain a single-stranded, positive-sense RNA ge¬ 
nome of »28 kb. Infection with these viruses causes wa¬ 
tery diarrhea, dehydration, and a high mortality rate among 
suckling pigs. Coronaviruses (CoVs) are prone to genetic 
evolution through accumulation of point mutations and ho¬ 
mologous recombination among members of the same ge¬ 
nus (7). Porcine respiratory coronavirus (PRCV), a mutant 
of TGEV, appeared in pigs in the 1980s (2). The spread 
of PRCV, which conserved most of the antigenic sites 
and causes cross-protection against TGEV (2), led to the 
gradual disappearance of TGEV. Newly emerging Co Vs 
pose a potential threat to human and animal health because 
multiple human CoV infections have been derived from 
animal hosts. Emerging swine coronaviruses are of great 
concern to swine health because of the potential increase 
in viral pathogenesis. 

In 1978, PEDV was first identified in Europe; subse¬ 
quent reports occurred in many countries in Asia, includ¬ 
ing China, Japan, Korea, and Thailand. In 2010-2012, 
genetically different PEDV strains emerged, causing se¬ 
vere outbreaks in China ( 4 ). PEDV spread to the United 
States, Canada, and Mexico in 2013-2014, and genetically 
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related strains were detected in South Korea and Taiwan 
(5-7). The PEDV outbreak caused large global economic 
losses to the swine industry. In Europe, a severe PEDV 
epidemic occurred in Italy during 2005-2006 (8), and 
in 2014-2015, PEDV was detected in Germany, France, 
and Belgium. These strains have a high nucleotide iden¬ 
tity to PEDV strains that contain distinct insertions and 
deletions (INDELs) in the S gene (S-INDELs) from the 
United States (9-11). We report the detection and genetic 
characterization of swine enteric Co Vs circulating in Italy 
during 2007-2014. We also report a recombinant TGEV 
and PEDV strain (identified as the species Swine enteric 
coronavirus [SeCoV]) circulating from June 2009 through 
2012. Finally, we describe the phylogenetic relationship of 
the 2014 PEDV S-INDELs to the recent PEDV strains cir¬ 
culating in Europe. 

The Study 

During 2007-2014, we collected 27 fecal and 24 intes¬ 
tinal samples from pigs with suspected PEDV or TGEV 
infections; the pigs came from swine farms in northern 
Italy (Po Valley), which contains the regions of Piemonte, 
Lombardia, Emilia Romagna, and Veneto (online Tech¬ 
nical Appendix Figure 1, http://wwwnc.cdc.gov/EID/ 
article/22/1/15-0544-Techappl.pdf). The Po Valley con¬ 
tains 70% of Italy’s swine. Clinical signs included watery 
diarrhea in sows and a death rate in piglets of 5%—10%, 
lower than is typical with PEDV or TGEV infections. Sam¬ 
ples were submitted for testing by electron microscopy, 
PEDV ELISA, viral isolation, pan-CoV reverse transcrip¬ 
tion PCR (RT-PCR), and RT-PCR for PRCV and TGEV; 
selected positive pan-CoV samples were sequenced (12- 
14) (online Technical Appendix). 

Results of electron microscopy showed that 25 (49%) 
of the 51 samples contained CoV-like particles, but all 
samples were negative for viral isolation. Although only 
38 samples (74%) were positive by pan-CoV RT-PCR, 
47 (92%) were positive by the PEDV ELISA (Table 1) 
(12,13). Of the 38 pan-CoV-positive samples, 18 were 
selected for partial RNA-dependent RNA polymerase 
(RdRp), spike (SI) (14), and membrane (M) sequencing 
(Table 1). All samples were negative for PRCV and TGEV 
by RT-PCR, ruling out co-infection with PEDV and TGEV 
or PRCV (15). 
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DISPATCHES 


Table 1. Distribution of test results of samples from pig farms in study of swine enteric coronaviruses in northern Italy, 2007-2014* 


Sample 

no. 

Farm 

no. 

Year 

Region 

EM 

PEDV 

ELISA 

Pan-CoV 

RT-PCR 

TGEV/ 
PRCV SI 

RdRp 

Sequences 

SI 

M 

222654 

1 

2007 

Emilia Romagna 

— 

+ 

+ 

NA 

NA 

NA 

NA 

1448 

2 

2007 

Emilia Romagna 

— 

+ 

— 

NA 

NA 

NA 

NA 

19908 

3 

2007 

Emilia Romagna 

— 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

70323 

4 

2007 

Lombardia 

+ 

+ 

+ 

NA 

NA 

t 

NA 

114372 

5 

2007 

Lombardia 

+ 

+ 

+ 

NA 

NA 

NA 

NA 

200079 

6 

2007 

Lombardia 

— 

+ 

+ 

NA 

NA 

t 

NA 

320855/5 

7 

2007 

Lombardia 

+ 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

320855/6 

7 

2007 

Lombardia 

+ 

+ 

+ 

NA 

t 

t 

NA 

3936/1 

8 

2008 

Lombardia 

— 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

3936/2 

8 

2008 

Lombardia 

— 

+ 

+ 

NA 

t 

NA 

NA 

29177 

9 

2008 

Veneto 

+ 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

43853 

10 

2008 

Lombardia 

+ 

+ 

— 

NA 

NA 

NA 

NA 

7239* 

11 

2009 

Lombardia 

— 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

20001 

12 

2009 

Lombardia 

— 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

20416 

13 

2009 

Lombardia 

— 

+ 

+ 

NA 

t 

t 

NA 

22603 

14 

2009 

Lombardia 

— 

+ 

+ 

— 

Cluster 1 

Cluster 1 

Cluster 1 

26199/2 

15 

2009 

Lombardia 

— 

+ 

— 

NA 

NA 

NA 

NA 

87565 

16 

2009 

Emilia Romagna 

— 

+ 

+ 

NA 

NA 

t 

NA 

111357/7 

17 

2009 

Lombardia 

NA 

+ 

— 

NA 

NA 

NA 

NA 

137442 

18 

2009 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

205396 

19 

2009 

Lombardia 

+ 

+ 

— 

NA 

NA 

NA 

NA 

208995 

20 

2009 

Lombardia 

+ 

— 

+ 

NA 

t 

NA 

NA 

213306* 

21 

2009 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

244945 

22 

2009 

Emilia Romagna 

+ 

+ 

— 

NA 

NA 

NA 

NA 

245242 

22 

2009 

Emilia Romagna 

+ 

+ 

+ 

NA 

t 

NA 

NA 

274771 

23 

2009 

Veneto 

+ 

— 

— 

NA 

NA 

NA 

NA 

307121 

24 

2009 

Emilia Romagna 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

315994 

25 

2009 

Lombardia 

+ 

— 

— 

NA 

NA 

NA 

NA 

320695 

26 

2009 

Lombardia 

+ 

+ 

+ 

NA 

NA 

t 

t 

320825 

26 

2009 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

324345 

27 

2009 

Lombardia 

+ 

+ 

+ 

NA 

t 

t 

t 

324374 

27 

2009 

Lombardia 

+ 

+ 

+ 

NA 

t 

t 

t 

324397 

27 

2009 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

324507/1 

28 

2010 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

324507/2 

28 

2010 

Lombardia 

+ 

+ 

+ 

NA 

NA 

t 

t 

324507/3 

28 

2010 

Lombardia 

+ 

+ 

+ 

NA 

NA 

NA 

t 

324507/4 

28 

2010 

Lombardia 

+ 

+ 

+ 

NA 

NA 

t 

t 

5448/2 

29 

2011 

Emilia Romagna 

NA 

+ 

— 

NA 

NA 

NA 

NA 

28607 

30 

2012 

Lombardia 

— 

+ 

+ 

NA 

NA 

t 

t 

29742 

30 

2012 

Lombardia 

+ 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

30917 

31 

2012 

Lombardia 

+ 

+ 

— 

NA 

NA 

NA 

NA 

35621/1 

32 

2012 

Lombardia 

+ 

+ 

— 

NA 

NA 

NA 

NA 

35621/2 

32 

2012 

Lombardia 

— 

+ 

+ 

NA 

NA 

NA 

NA 

41906 

33 

2012 

Lombardia 

— 

+ 

+ 

NA 

NA 

NA 

NA 

44833 

34 

2012 

Lombardia 

NA 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

67322 

8 

2012 

Lombardia 

— 

+ 

+ 

NA 

NA 

NA 

t 

273992 

35 

2012 

Lombardia 

— 

+ 

+ 

— 

Cluster II 

Cluster II 

Cluster II 

32961 

36 

2013 

Piemonte 

— 

+ 

— 

NA 

NA 

NA 

NA 

32963 

36 

2013 

Piemonte 

+ 

+ 

— 

NA 

NA 

NA 

NA 

178509 

37 

2014 

Emilia Romagna 

NA 

NA 

+ 

— 

Cluster III 

Cluster III 

Cluster III 

200885 

38 

2014 

Emilia Romagna 

+ 

+ 

+ 

— 

Cluster III 

Cluster III 

Cluster III 


*Cluster I represents strains circulating from 2007 through mid-2009; cluster II represents strains circulating from mid-2009 through 2012; and cluster III 
represents strains circulating since 2014. EM, electron microscopy; M, membrane; pan-CoV RT-PCR, pan-coronavirus reverse transcription PCR; PEDV, 
porcine epidemic diarrhea virus; PRCV, porcine respiratory coronavirus; RdRp, RNA-dependent RNA polymerase; SI, spike; TGEV, transmissible 
gastroenteritis virus; +, positive; -, negative; NA, not tested or sequenced. 
tSequence available but not included in this study. 

^Samples selected for whole genome sequencing._ 


On the basis of the partial sequences from RdRp and 
the S1 and M genes, the strains from Italy clustered into 3 
temporally divided groups, suggesting 3 independent virus 
entries. Cluster I represents strains circulating from 2007 
through mid-2009; cluster II represents strains circulating 
from mid-2009 through 2012; and cluster III represents 


strains circulating since 2014 (online Technical Appendix 
Figure 2, panels A-C). Cluster I was identified in Emilia Ro¬ 
magna (n = 1), Lombardia (n = 5), and Veneto (n = 1). Clus¬ 
ter II was identified in Emilia Romagna (n = 1) and Lom¬ 
bardia (n = 8). Cluster III was identified in Emilia Romagna 
at 2 swine farms. To help explain the temporal clustering, a 
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single S1 gene segment was sequenced from clusters I and 
II (PEDV/Italy/7239/2009 and SeCoV/Italy/213306/2009, 
respectively). Because of the recent outbreak of PEDV in 
Europe, the 2 positive samples from cluster III (PEDV/ 
Italy/178509/2014 and PEDV/Italy/200885/2014) were se¬ 
quenced (Figure 1, panel A). 

One strain from each cluster was selected for whole 
genome sequencing (online Technical Appendix). Un¬ 
fortunately, the whole genome was obtained from only 
clusters I and II (PEDV/Italy/7239/2009 and SeCoV/Ita- 
ly/213306/2009, respectively; Figure 1, panel B). Recombi¬ 
nation analysis was conducted on the 2 whole genomes and 
was not detected in PEDV/Italy/7239/2009. Recombina¬ 
tion was detected in SeCoV/Italy/213306/2009 at position 
20636 and 24867 of PEDV CV777 and at position 20366 
and 24996 of TGEV HI6 (Figure 2), suggesting the oc¬ 
currence of a recombination event between a PEDV and a 
TGEV. The complete S gene of SeCoV/Italy/213306/2009 


shared 92% and 90% nt identity with the prototype Euro¬ 
pean strain PEDV CV777 and the original highly virulent 
North American strain Colorado 2013, respectively, and 
the remaining genome shared a 97% nt identity with the 
virulent strains TGEV HI6 and TGEV Miller M6. Whole- 
genome analysis of PEDV/Italy/7239/2009 showed that it 
grouped with the global PEDV strains (<5) and shared ~97% 
nt identity with PEDV strains CV777, DR 13 virulent, and 
North American S-INDEL strain OH851 (Table 2). 

Conclusions 

During 2007-2014, most (92%) samples collected from 
the Po Valley in Italy were positive for PEDV by ELI¬ 
SA; only 72% were positive by pan-CoV PCR. However, 
because we were investigating the presence of PEDV 
or TGEV in samples with clinical signs of diarrhea, the 
high occurrence of PEDV may not reflect the actual 
prevalence of PEDV in Italy. The increased percentage 
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Figure 1. Phylogenetic analyses of swine enteric coronaviruses in Italy. A) Analysis performed on the basis of the nucleotide sequence 
of the complete spike (SI) gene of 4 representative strains from the 3 clusters and B) whole genome of 2 positive strains from clusters 
I and II. Cluster I represents strains circulating from 2007 through mid-2009; cluster II represents strains circulating from mid-2009 
through 2012; and cluster III represents strains circulating since 2014. Bootstrap values >70% (1,000 replicates) are indicated. 
Reference sequences are identified by GenBank accession no. and strain name. The strains from this study are represented in blue; 
strains from China are green; strains from North America are red; and strains from Europe are purple. PEDV, porcine epidemic diarrhea 
virus; PRCV, porcine respiratory coronavirus; TGEV, transmissible gastroenteritis virus; SeCoV, swine enteric coronavirus. 
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DISPATCHES 
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Figure 2. Potential recombination points in the SeCoV strains in study of swine enteric coronaviruses in Italy. The potential parent 
strains HI6 (TGEV) and CV777 (PEDV) are shown in teal and purple, respectively. Arrows indicate recombinant breakpoints. UTR, 
untranslated region; ORF, openreading frame; S, spike; E, envelope; M, membrane; N, nucleocapsid; nsp, nonstructural protein; PEDV, 
porcine epidemic diarrhea virus; TGEV, transmissible gastroenteritis virus; SeCoV, swine enteric coronavirus. 


of PEDV found in samples tested by ELISA, compared 
with the proportion found by PCR, may be explained by 
the number of ambiguous bases in the pan-CoV primers; 
the ambiguous bases severely reduce the efficiency of the 
reaction. The swine enteric CoV strains from Italy in our 
study, including the recombinant strain, were reported in 
pigs with mild clinical signs, indicating that PEDV and 
SeCoV have been circulating in Italy and likely through¬ 
out Europe for multiple years but were underestimated as 
a mild form of diarrhea. 

To understand the evolution of PEDV in Italy, the par¬ 
tial RdRp, S, and M genes were sequenced from 18 sam¬ 
ples and grouped in 3 different temporal clusters. Cluster I 


(2007-mid 2009) resembles the oldest PEDV strains; clus¬ 
ter II resembles a new TGEV and PEDV recombinant vari¬ 
ant; and cluster III, identified from 2 pig farms in northern 
Italy in 2014, resembles the PEDV S-INDEL strains identi¬ 
fied in Germany, France, Belgium, and the United States. 
The >99.3% nt identity of the SI gene within cluster III and 
in previously identified strains could suggest the spread of 
the S-INDEL strain into Europe. However, directionality 
of spread cannot be determined because of a lack of global 
and temporal PEDV sequences. 

Although our findings could indicate 3 introductions 
of PEDV in Italy, the results more likely suggest the high 
ability of natural recombination among Co Vs and the 


Table 2. Nucleotide identities of strains PEDV/ltaly/7239/2009 and SeCoV/ltaly/213306/2009, representative of clusters I and II, 
respectively, in study of swine enteric coronaviruses in Italy* _ 

ORF1 Spike ORF3 Envelope Membrane Nucleocapsid 


Strain identification 

1 

II 

1 

II 

1 

II 

1 

II 

1 

II 

1 

II 

PEDV/Belgium/CV777/1977 

97.3 

57.8 

96.3 

92.7 

98.1 

43.6 

97.6 

43.0 

97.6 

55.9 

96.6 

42.5 

PEDV/South Korea/DR13 vir/2009 

98.1 

58.0 

97.3 

93.1 

99.3 

43.7 

99.6 

43.0 

97.7 

55.6 

97.6 

43.3 

PEDV/USA/Colorado/2013 

98.0 

57.9 

94.6 

90.7 

98.6 

44.0 

99.2 

43.0 

97.8 

55.3 

96.8 

43.4 

PEDV/USA/OH851/2014 

98.1 

57.9 

96.9 

91.9 

98.7 

43.8 

99.2 

43.0 

97.7 

55.5 

96.7 

43.4 

PEDV/L00721/GER/2014 

98.0 

57.9 

97.0 

92.0 

98.6 

43.7 

99.2 

43.0 

97.8 

55.3 

96.8 

43.4 

TGEV/USA/Miller M6/2006 

57.9 

96.8 

52.0 

52.5 

40.0 

89.1 

42.6 

96.4 

56.1 

97.1 

43.1 

96.3 

PRCV/USA/ISU-1/2006 

58.0 

96.5 

47.7 

48.2 

53.0 

76.6 

43.8 

96.0 

56.2 

96.3 

42.8 

95.7 


*Cluster I represents strains circulating from 2007 through mid-2009; cluster II represents strains circulating from mid-2009 through 2012. Ger, Germany; 
ORF, open reading frame; PRCV, porcine respiratory coronavirus; PEDV, porcine epidemic diarrhea virus; SeCoV, swine enteric coronavirus; TGEV, 
transmissible gastroenteritis virus. 
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continued emergence of novel Co Vs with distinct patho¬ 
genic properties. Further investigation is needed to deter¬ 
mine the ancestor of the SeCoV strain or to verify whether 
the recombinant virus was introduced in Italy. Recombinant 
SeCoV was probably generated in a country in which both 
PEDV and TGEV are endemic, but because the presence of 
these viruses in Europe is unclear and SeCoV has not been 
previously described, it is difficult to determine the paren¬ 
tal strains and geographic spread of SeCoV. Future studies 
are required to describe the pathogenesis of SeCoV and its 
prevalence in other countries. 
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